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I. INTRODUCTION 



The striking features of the nucleon structure func- 
tion i<2 were first noted nearly 30 years ago by Bloom 
and Gilman pj. They empirically observed two effects 
in data measured at SLAC: a) the dual behaviour of 
the F2(x,Q 2 ) function that shows common features be- 
tween the two kinematic regions corresponding to the 
nucleon resonances and Deep Inelastic Scattering (DIS), 
b) the extension of the scaling region to lower Q 2 val- 
ues when F2(x',Q 2 ) is plotted as a function of x' — 
x/(l + M 2 x 2 /Q 2 ),the "improved scaling variable" . More 
precisely, they found that the smooth function F 2 (x') 
measured at high Q 2 in the DIS region represents a good 
average over the resonances of the F2(x',Q 2 ) structure 
function measured at lower Q 2 values. Moreover, the 
duality appears to be valid locally. In fact, each of the 
most prominent resonance bumps, when averaged within 
its width, shows approximate scaling Q • Later on, in the 
framework of QCD, De Rujula, Georgi and Politzer Q 
provided the first explanation of the Bloom-Gilman du- 
ality. They evaluated the Cornwall-Norton Q moments 
(Q 2 ) of the nucleon structure function F2, defined 
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as 

M^ N (Q 2 )= f 1 dxx (n - 2 ^F 2 (x,Q 2 ), (1) 
Jo 

using the Operator Product Expansion (OPE) they ob- 
tained the following expression 

00 / 2 \ & 

M% N (Q 2 ) = A n (Q 2 ) + J2y^) B nk (Q 2 ), (2) 

where A n (Q 2 ) can be evaluated in the framework of per- 
turbative QCD (pQCD), and it is directly connected to 
the corresponding moment of the asymptotic limit of F 2 . 
The contribution of the higher twists, which is related to 
multi-parton correlations inside the nucleon and repre- 
sented by B n k(Q 2 ), depends on the value of the constant 
7 2 in such a way that j 2 can be considered as a scale 
constant for higher twist effects. Assuming a small value 
of the constant 7 2 , the authors of Ref. Q showed that 
the contribution of the higher twists was relatively small, 
at least for low values of n and for Q 2 > M 2 , justifying 
the observed dual behaviour of the structure function. 

It is now well established that the interpretation of 
the parton-hadron duality in light of QCD requires the 
evaluation of the moments of the nucleon structure func- 
tions and their evolution as a function of Q 2 . Current 
pQCD calculations can estimate the Q 2 evolution up to 
the Next-to-Next-to-Leading Order, giving access to the 
interesting kinematic region of high x and moderate Q 2 
where the multi-parton correlation contribution to the 
nucleon wave function becomes dominant. 
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FIG. 1: Experimental data on the structure function 
F2{x,Q 2 ) used for the moment evaluation in the CLAS kine- 
matic region: points - world data; shaded area - CLAS data. 



The interest in investigating multi-parton correlations 
in inelastic lepton scattering off the nucleon at large val- 
ues of x has recently been renewed, leading to a re- 
analysis of old F 2 data 0, |(| . Unfortunately, the results 
from Ref. and were mainly based on the analysis 
of fits of the structure function F 2 and therefore were 
still qualitative. Moreover, the previous lack of data in 
the resonance region did not allow a model independent 
evaluation of the moments' evolution to lower Q 2 , and 
therefore offered very few opportunities to quantitatively 
investigate the role of QCD below the DIS limit. 

The Hall C Collaboration at the Thomas Jefferson Na- 
tional Accelerator Facility (TJNAF) has recently pro- 
vided high quality data in this kinematic region , allow- 
ing a more precise evaluation of the moments of the F 2 
structure function of the proton Q. However, like many 
other such measurements, the data were taken with a 
spectrometer of relatively small angular acceptance and 
the measured inclusive cross sections do not span a large 
continuous x interval for constant Q 2 . Data taken in this 
manner follow a kinematic locus in Q 2 vs. x and require 
substantial interpolation to determine the F 2 moments. 

In this paper we report the first measurement in a wide 
continuous interval in x and Q 2 (see Fig. of the inclu- 
sive electron-proton scattering cross section. These mea- 
surements were performed at TJNAF with the CLAS de- 
tector in Hall B. The F 2 structure function was extracted 
over the whole resonance region (W < 2.5 GeV) below 
Q 2 = 4.5 (GeV/c) 2 . This measurement, together with 
existing world data, allowed for the evaluation of the F 2 
moments, drastically reducing the uncertainties related 
to data interpolation and providing the most detailed 
dependence on Q 2 of the moments up to n = 8. Fur- 
thermore, the elastic contribution to the moments was 
updated with respect to Ref. 13 using the fit of the nu- 
cleon form factors from Ref. |8| adjusted to the Jeffer- 



son Lab data on the ratio Ge/Gm 3, as described in 
Ref. |ldj| . Finally, we used our new determination of the 
F 2 moments to extract the higher twist contribution as 
a function of Q 2 . 

In section HT1 we review the F 2 moments in the frame- 
work of pQCD. In section IIHI we discuss the analysis 
of the data, including the extraction of the F 2 struc- 
ture function from the cross section. The evaluation of 
the moments and uncertainties is also presented in sec- 
tion ^Q] Finally, Section Hvl is devoted to the interpreta- 
tion of the results. 



II. MOMENTS OF THE STRUCTURE 
FUNCTION F 2 

Until recently the studies of inclusive lepton-nucleon 
scattering represented the main source of information 
about nucleon structure. In the DIS region, measured 
structure functions can be directly connected to the par- 
ton momentum distribution of the nucleon in the frame- 
work of pQCD. After the successful interpretation of the 
DIS region, the intermediate kinematic domain, situated 
at Q 2 of a few (GeV/c) 2 and large values of x, attracted 
the interest of physicists. Despite interpretation difficul- 
ties, this region allows the study of multi-parton corre- 
lation contribution to the proton wave function. These 
processes are not accessible in DIS due to the small value 
of the running coupling constant as{Q 2 )- 

The OPE of the virtual photon-nucleon scattering am- 
plitude leads to the description of the complete Q 2 evolu- 
tion of the moments of the nucleon structure functions. 
The n-th Cornwall-Norton moment [j] of the (asymp- 
totic) structure function F 2 {x, Q 2 ) for a massless nucleon 
can be expanded as: 



M n (Q ) 



T=2k 



§0-2) 



(3) 



where fc = l,2,...,oo,/iis the factorization scale, O nT (fi) 
is the reduced matrix element of the local operators with 
definite spin n and twist r (dimension minus spin), re- 
lated to the non-perturbative structure of the target. 
EnrilJ'i Q 2 ) is a dimensionless coefficient function describ- 
ing the small distance behaviour, which can be pertur- 
batively expressed as a power expansion of the running 
coupling constant a s (Q 2 ). 

At Q 2 values comparable with the squared proton 
mass, M 2 , the structure function F 2 still contains non- 
negligible mass-dependent terms that produce in Eq. [3] 
additional M 2 /Q 2 power corrections (kinematic twists). 
To avoid these terms, the moments M^ N {Q 2 ) of the 
massless F 2 have to be replaced in Eq. [3] by the corre- 
sponding Nachtmann moments M^(Q 2 ) of the mea- 
sured structure function F 2 {x,Q 2 ) (see also Ref. ^ 
has been shown that: 



MZ N (F» m (x,Q 2 )) 



M»(F 2 (x,Q 2 )), 



(4) 
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III. DATA ANALYSIS 





FIG. 2: Twist diagrams: (a) the leading twist contribution 
evaluated at leading order of pQCD; (b) the contribution of 
higher twists, where current quark and nucleon remnant can 
exchange by a system of particles consisting of gluons and 
qg-pairs whose complexity is increasing with twist order. 



where ^""(x, Q 2 ) is the asymptotic structure function 
of the massless nucleon and 



cn+l 

dx^ r F 2 {x,Q 2 ) 
x 3 

3 + 3(n + l)r + n(n + 2)r 2 
(n + 2)(n + 3) 



(5) 



where r = y/l + AM 2 x 2 /Q 2 and £ = 2i/(l + r). 

Since the moments in Eq. [3] are totally inclusive, the 
elastic contribution at x = 1 has to be added according 
toRef. 0: 



F% l (x,Q 2 ) = 5(l 



x) 



(G|(0 2 



■jjrzG 2 M (Q 2 )) 



1 



(0) 



liii- 



with G 2 E {G 2 M ) being the proton electric (magnetic) elas- 
tic form factor. 

For the leading r = 2 twist, one ends up in Leading 
Order (LO) or Next-to-Leading Order (NLO) with the 
well known perturbative logarithmic Q 2 evolution of sin- 
glet and non-singlet F2 moments. However, if one wants 
to extend the analysis to small Q 2 and large x where 
the rest of the perturbative series becomes significant, 
some procedure for the summation of the higher orders 
of the pQC D expansion, such as the infrared renormalon 
model ]l3LIT^ | or the recently developed soft-gluon resum- 
mation technique |l5l Il6|. has to be applied. For higher 
twists, r > 2, the power terms E nT (p,) are related to 
quark-quark and quark-gluon correlations, as illustrated 
by Fig. [3 and should become important at small Q 2 . 

The systematic analysis of the Q 2 dependence of the 
experimentally derived Nachtmann moments (Q 2 ) in 
the intermediate Q 2 range (0.5 < Q 2 < 10 (GeV/c) 2 ) 
should allow a separation of the higher twists from the 
leading twist. A precise evaluation would permit a com- 
parison with the QCD predictions obtained from lattice 
simulations or a comparison with those models that de- 
scribe the non-perturbative domain. 



The data were collected at TJNAF in Hall B with the 
CLAS detector and a liquid hydrogen target with thick- 
ness px = 0.35 g/cm 2 during the electron beam running 
period in February-March 1999. The average beam cur- 
rent of 4.5 nA corresponded to a luminosity of 6xl0 33 
cm -2 s . To cover the largest interval in Q 2 and x, 
data were taken at five different electron beam energies: 
E = 1.5, 2.5, 4.0, 4.2 and 4.4 GeV. The CLAS detec- 
tor is a magnetic spectrometer based on a six-coil torus 
magnet whose field is primarily oriented along the az- 
imuthal direction. The sectors, located between the mag- 
net coils, are instrumented individually to form six es- 
sentially independent magnetic spectrometers. The par- 
ticle detection system includes drift chambers (DC) for 
track reconstruction 01 > scintillation counters (TOF) for 
the time of flight measurement B3, Cherenkov counters 
(CC) for electron identification [n|, and electromagnetic 
calorimeters (EC) to measure neutrals and to improve 
the electron-pion separation 20] . The EC detectors have 
a granularity defined by triangular cells in the plane per- 
pendicular to the incoming particles to study the elec- 
tromagnetic shower shape and are longitudinally divided 
into two parts with the inner part acting as a pre-shower. 
Charged particles can be detected and identified for mo- 
menta down to 0.2 GeV/c and for polar angles between 
8° and 142°. The CLAS superconducting coils reduce 
the acceptance of about 80% at 9 = 90° to about 50% 
at forward angles (9 — 20°), while the total acceptance 
for electrons is about 1.5 sr. Electron momentum reso- 
lution is a function of the scattered electron angle and 
it varies from 0.5% for 9 < 30° up to 1-2% for 9 > 30°. 
The angular resolution is approximately constant and ap- 
proaching 1 mrad for polar and 4 mrad for azimuthal 
angles: the resolution on the momentum transfer ranges 
therefore from 0.2 up to 0.5 %. The missing mass resolu- 
tion was estimated 2.5 McV for beam energy less than 3 
GeV and about 7 MeV for larger energies. To study all 
possible multi-particle production, the acquisition trig- 
ger was configured to require at least one electron can- 
didate in any of the sectors, where an electron candidate 
was defined as the coincidence of a signal in the EC and 
Cherenkov modules for each sector separately. 

The accumulated statistics at the five energies is large 
enough (> 6 x 10 s triggers) to allow the extraction of 
the inclusive cross section with a rather small statisti- 
cal error (< 5%), in small x and Q 2 bins (Ax = 0.009, 
AQ 2 = 0.05 GeV 2 ). The determination of the systematic 
error was more critical. CLAS is a large acceptance spec- 
trometer and the response depends on the energy E' and 
the angle 9 of the scattered electron. Determining the 
systematic effects of these, and other experimental pa- 
rameters, is both necessary and complex. Consequently, 
we dedicate the next subsections to the discussion of the 
analysis procedure. 



5 



A. Momentum Correction 

Determining the momentum of a charged particle mea- 
sured with CLAS depends on a proper understanding of 
the magnetic field geometry. Due to the complexity of 
the detector and particularly the torus magnet system, 
it is crucial to check the reliability of the momentum de- 
termined by the DC tracking system. For this reason the 
position of the elastic peak was extracted from the mea- 
sured inclusive electron cross section and compared to 
the theoretical value. A correction to the scattered elec- 
tron momentum was applied to shift the elastic peak to 
the accepted value. The momentum correction obtained 
was small (from 2 to 7 MeV in W, depending on the 
beam energy) and resulted in significant improvement in 
the width of elastic peak . 

The systematic error on the correction was estimated 
by comparing the position of the well known second res- 
onance_(Sjj(1535) + Di3(1520)) to the position given in 
Ref. Hit |22j. The position difference AW affects the 
cross section evaluation. The relative systematic error 
on the momentum correction is therefore given by 



\a B {W 



*f,Q*)-a B (W+^f,Q 2 ) 



<r B (W,Q 2 



(7) 



where a B represents the Bodek fit according to the 
parametrization from Ref. 0, ■ The systematic er- 
ror calculated with Eq. is given in Table. D 



B. Electron Identification and Pion Rejection 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
ECtot/P 

FIG. 3: Photoelectron distribution in the Cherenkov detector 
versus the energy deposited in the EC detector divided by the 
momentum of the particle as determined by the drift chamber. 
The black vertical line represents the cut to reduce the pion 
contamination. 



and only the pion contamination to the left of the solid 
line can be removed without affecting the electron detec- 
tion efficiency. The remaining pion contamination and 
the correction of the Cherenkov efficiency for electrons 
Fphe have determined by a combined fit of the measured 
photoelectron distribution with two Poisson distributions 
convoluted with a Gaussian function to account for the 
finite photomultiplier resolution as shown in Fig. 0] 



The electrons were identified by a combined off-line 
analysis of the signals from the four detector systems 
(DC, TOF, CC and EC). Only those electron candi- 
dates that were detected inside the most uniform ( "fidu- 
cial" ) detector volume were analyzed. The electron yield 
was corrected separately in each kinematic bin, for pion 
contamination, detection efficiency and radiative correc- 
tions. 

The photoelectron distribution in the CC depends on 
the kinematics, and the contaminating pion peak can be 
completely removed only with large efficiency losses of 
about 30% in several kinematic regions. Therefore the 
pion contamination was removed by a two-step proce- 
dure. Electrons producing a large number of photoelec- 
trons (see Fig. were identified by an energy cut in the 
EC detector response. The pion contamination to elec- 
trons producing a small number of photoelectrons was 
removed by analysis of the CC response. 

As an example of the first step, Fig. [3| shows the CC 
photoelectron distribution N p h e as a function of the frac- 
tion of energy deposited in the EC detector EC to t/P 
for negatively charged particles emitted at 9 < 35°, mo- 
mentum P < 1 GeV/c, and a beam energy of 2.5 GeV. 
The regions corresponding to pions (N p h e < 2.5) and to 
electrons (ECtot/P > 0.25) cannot be clearly separated 
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FIG. 4: The fitted photoelectron distribution for two different 
sets of kinematics after removing some of the pion contamina- 
tion via the EC tot / P cut: solid triangles show the distribution 
obtained with 4 GeV beam, scattered electron angle 9 = 31° 
and momentum P = 1 GeV; empty diamonds represent data 
taken with 2.5 GeV beam, scattered electron angle 9 — 41° 
and momentum P = 1 GeV. 
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The fit was performed separately for each sector over 
the whole kinematics data set (8 = 20° — 50° and W — 
0.9 - 2.5 GeV) 

To minimize the errors, the fit was performed in rather 
large bins (A9 = 2° and AE 1 ~ 0.1 GeV). Therefore, 
in order to apply the correction to the measured cross 
section, which was obtained with smaller bins, values of 
the correction were parametrized with the polynomial 
function: 

F phe (v,0) =l + A(y- v ) + {B + C6)(y - i> ) 2 , (8) 

where A, B, C and i>q are free parameters and v = Eq—E' 
the electron energy transfer. The related systematic error 
Sphe — SF p he/ F p he is mainly due to the low statistics 
in those bins corresponding to large Q 2 values and was 
found to be less than 2%. 



in detail including all materials that make up each detec- 
tor. The simulations of detector response to the scattered 
electron were performed according to the following pro- 
cedure: 

1. Electron scattering events were generated by a 
random event generator with the probability dis- 
tributed according to <rf ad , described above. The 
values of elastic and inelastic cross sections of the 
electron-proton scattering were taken from exist- 
ing fits of world data, in references Q and 0, , 
respectively. The internal radiative processes con- 
tribution was added according to calculations [25| . 

2. The generated events were passed through the 
standard CLAS GEANT-based simulation pro- 
gram 26], to model the detector response. 



C. Background Subtraction 

Since the pair production background has not been 
measu red, its contribution was estimated according to a 
model ,23 based on the Wiser fit [24( of the inclusive pion 
photo-production reaction. The most important source 
of e + e~ pairs in the CLAS is due to ir° production, which 
either decays to -fe + e~~ (Dalitz decay) or to 77, with 
subsequent photon conversion to e + e~. The model was 
carefully checked, and it was in good agreement with 
the measured positron cross sections |23j; the difference 
was always less than 30%. The value of the correction 
was assumed to be equal to the ratio of the inclusive 
e + production cross section a e + over Bodek's fit 0, 
including radiative processes (tail from the elastic peak, 
bremsstrahlung, and Schwinger correction) crf ad : 



F e + e - (Eq, E', 9) — 



1 



cr e + (Eq.E' , 
a? ad (E ,E> 



(9) 



The correction is generally small as expected in Ref. 7] , 
therefore it was applied only for E n > 2.0 GeV and W > 
1.7 GeV, where it was about 2%. The relative systematic 
error 5 e + e - = SF e + e -/F e + e - from this correction was 
estimated using uncertainties on cr e + given in Ref. |23j . 

In order to remove the contribution of scattering on 
the target walls, the empty target data were analyzed 
in the same way and subtracted from the inclusive data, 
after proper normalization. An additional source of back- 
ground originating from knock-on electrons produced in 
the supporting structure of the detector was estimated 
and it found to be smaller than 0.3%. 



3. The results of the previous stage were further pro- 
cessed to make the detector response more realis- 
tic by adding the effects of electronic noise, back- 
ground, dead wires and scintillator paddles. 

4. Finally, the efficiency was calculated in each kine- 
matic bin as a ratio of the number of reconstructed 
events N rec over the number of generated events 



N, 



gen ■ 



e e ff{E ,x,Q 2 



N rec (E ,x,Q 2 ) 
N gen (E ,x,Q 2 ) 



(10) 



The electron detection efficiency obtained from simula- 
tions is about 97% and approximately constant inside 
the fiducial region of the detector over the whole avail- 
able kinematics. 

In order to test the reliability of the simulation proce- 
dure and to check a proper absolute normalization, the 
well known elastic scattering cross section was extracted 
from the same data set (^ exp ) and compared to the 

simulated cross section (4ft . ). The two cross sections 
are in good agreement within statistical and systematic 
errors as shown in Fig. |3J 

The relative systematic deviation of the elastic cross 
section obtained from simulations and from these data 
S e ff, was calculated for each beam energy and scattered 
electron angle (in bins of one degree on the accessible 
interval from 20 to 50 degrees) according to: 



S 2 e ff (Eo,0) 



+ S 2 exp (E ,6) = 



da I rp 



(11) 



D. Simulations 

Due to the complexity of the CLAS detector the only 
way to study its response functions is to perform com- 
plete computer simulations, describing each subsystem 



where ^j it is the parametrization described in Ref. @ 
and S exp is the statistical error of the measured elas- 
tic cross section. For the error propagation S e ff was 
parametrized by a linear function of the scattered elec- 
tron angle 9. 
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FIG. 5: Typical ratio of the measured elastic scattering cross 
section to the parametrization from Refs. 0.0] (points) with 
radiative corrections, in comparison to that obtained from the 
simulations (continuous line); errors are statistical only. 



All systematic uncertainties were propagated in 
quadrature to the final relative systematic error: 

S sys (x, Q 2 ) = [S 2 eff (x, Q 2 ) + S 2 phe (x, Q 2 ) + 

S 2 e+e - (x, Q 2 ) + S 2 mom (x, Q 2 ) + 5 2 rad (x, Q 2 )} 1/2 ,(15) 

where 5 ra d is the systematic uncertainty on the radiative 
correction, given by: 



Srad(E, x, Q 2 ) 



,TSAI 



(E,x, 



POLRAD 



(E,x,Q 2 ) 



where ei'^f- 1 (E , x , Q 2 



and 



.POLRAD 



(16) 

(E, x, Q 2 ) are the 



TSAI I 

e rad l^^iV J c rad 

radiative correction factors in cr^ ad evaluated with two 

[27|). These two ^approaches 
" and EF 



different approaches (|2E 

use different parametrizations of the elastic and |22j) 
and inelastic ([2ll 1221] and [2!|) cross sections as well as 
different calculation techniques. 

S r ad(E, x, Q 2 ) varies in the kinematic range of the ex- 
periment from to 20% while the average value is 3%. 
A minimum radiative correction systematic error of 1.5% 
was assumed. 



E. Inclusive Inelastic Cross Section 



F. Structure Function F2(x,Q 2 ) 



Since the Monte Carlo simulations were shown to be re- 
liable, they were used to evaluate efficiency, acceptance, 
bin centering and radiative corrections. For each kine- 
matic bin, the inclusive cross section da and the structure 
function F 2 were extracted directly from the raw electron 
yield N exp normalized to the integrated luminosity and 
corrected for efficiency, acceptance, bin centering, and 
radiative effects as follows: 



d 2 a 
dxdQ 2 



1 



N exp {x,Q 2 ) 



P^EQtot <x,Q 2 ) 



(12) 



where p is the density of liquid H 2 in the target, Na is 
the Avogadro constant, Ma is the target molar mass, L is 
the target length, Q to t is the total charge in the Faraday 
Cup (FC) and e(x, Q 2 ) is the efficiency defined in Eg. 1101 
with the radiative and bin-centering correction factors 
according to: 



e(x, Q 2 ) = e e ff{x, Q 2 )e rad {x, Q 2 )e bin {x, Q 2 ) , (13) 



where 



^rad 



and Ebi 



(14) 



and the integral was taken over the current bin area At. 
The radiative correction factor e ra d strongly varies in the 
explored kinematic range from 0.85 up to 1.6. Fortu- 
nately, the largest correction was contributed by the elas- 
tic peak tail for which calculations are very accurate (see 
Refs. jHlil). 



The structure function F 2 (x,Q 2 ) was extracted from 
the inelastic cross section using the fit of the function 
R(x, Q 2 ) = <jl/(Tt developed in [lj] and described in Ap- 
pendix E] The inclusive electron scattering cross section 
can be expressed in terms of the well known structure 
functions W\ and W 2 as [T^ : 



- =a Mo «{21F 1 (x,Q 2 )tan 2 - + VF 2 (x,Q 2 )} , (17) 



dD.dE' 1 x ' 2 

where the Mott cross section is given by: 

a 2 cos 2 % 

&Mott 



AE 2 sin 4 



(18) 



Therefore, the structure function F% — vWi is given by: 



F 2 (x,Q 2 ) = 



1 d 2 a 



a Mott dxdQ 2 l + l=£_i_ ' 



where J is the Jacobian given by 

_ xjs - M 2 ) 
~ 2-kMv ' 



(19) 



(20) 



where s is the squared invariant mass of the initial 
electron-proton system s = M 2 + 2EM and e is the po- 
larization parameter defined as 



e= l + 2 



v 2 + Q\ .. i 



2 tan 2 



(21) 



The function R(x,Q 2 ) is poorly known in the reso- 
nance region, however the structure function F 2 in the 



8 



relevant kinematic range is very insensitive to the value 
of R. In fact even a 100% systematic uncertainty on R 
gives only a few percent uncertainty on F 2 . The relative 
total systematic error is given by: 



6t s (x,Q 2 ) + 



1/2 



heRl + RI 

(22) 

The uncertainties of R given in Ref. Q were propagated 
to the resulting F 2 , and the actual systematic errors in- 
troduced by 5r were always lower than 3%. 

The combined statistical and systematic precision of 
the obtained structure function F 2 is strongly depen- 
dent on kinematics and the statistical errors vary from 
0.2% up to 30% at the largest Q 2 where statistics are 
very limited. Fig. shows a comparison between the 
F2 data from CLAS and the other world data in the 
Q 2 = 0.775 GeV 2 bin. The observed discrepancies with 
the data from Ref. Q which fill the large x region in 
Fig. are mostly within the systematic errors. Because 
of the much smaller bin centering corrections in this Q 2 
region our data are in a better agreement with data pre- 
viously measured at SLAC, given in Ref. , and the 
parameterization of those from Ref. [101 E2- The average 
statistical uncertainty is about 5%; the systematic uncer- 
tainties range from 2.5% up to 30%, with the mean value 
estimated as 7.7% (see Table|TJ). The values of F 2 (x,Q 2 ) 
determined using our data are tabulated elsewhere |lCj . 



TABLE I: Range and average of systematic errors on F2 



Source of uncertainties 


Variation range 


Average 


l%\ 


\%\ 


Efficiency evaluation 


1-9 


4.3 


e + e~ pair production correction 


0-3 


0.3 


Photoelectron correction 


0.1-2.2 


0.6 


Radiative correction 


1.5-20 


3.2 


Momentum correction 


0.1-30 


3.5 


Uncertainty of R = ^Jj. 


0.5-5 


2.4 


Total 


2.5-30 


7.7 



G. Moments of the Structure Function F-z 

As discussed in the introduction, the final goal of this 
analysis is the evaluation of the Nachtmann moments of 
the structure function F 2 . The total Nachtmann mo- 
ments were computed as the sum of the elastic and in- 
elastic moments: 



Ml 



(23) 



The contribution originating from the elastic peak was 
calculated according to the following expression from 




FIG. 6: Structure function F 2 (x,Q 2 ) at Q 2 = 0.775 GeV 2 : 
stars represent experimental data obtained in the present 
analysis with systematic errors indicated by the hatched 
area, empty circles show data from previous experiments 

BSIalEilMHIMISISHIM and 



Ref. \U 



Mi 



n+l 



3 + 3(n + l)r + n(n + 2)r 2 
(n + 2)(n + 3) 



GUQ 2 ) + S^G 2 M (Q 2 ) 



1 



Q 2 

4M- 



(24) 



where the proton form factors G%(Q 2 ) and G 2 M (Q 2 ) are 
from Ref. |8( modified according the rec ently measured 
data on Ge/Gm Q as described in Ref. |10|. 

The evaluation of the inelastic moment M" 1 involves 
the computation at fixed Q 2 of an integral over x. For 
this purpose, in addition to the results obtained from the 
CLAS data, world data on the structure function F 2 from 

Refsji m m in m m m m m m m m m pi 

I43ll44j and data on the inelastic cross section |2lll22Ll45| 

were used to reach an adequate coverage (see Fig. 0. 
The integral over x was performed numerically using the 
standard trapezoidal method TRAPER ^(|. Data from 
Ref. ^3 were not included in the analysis due to their 
inconsistency with other data sets as explained in detail 
in Ref. |48|. and data from Ref. 0,H(j were not included 
due to the large experimental uncertainties. 

The Q 2 -range from 0.05 to 3.75 (GeV/c) 2 was divided 
into AQ 2 = 0.05 (GeV/c) 2 bins. Then within each Q 2 
bin the world data were shifted to the central bin value 
Ql, using the fit of F^(x, Q 2 ) from Ref. Q. Here the fit 
Fn *(x, Q 2 ) consists of two parts, a parametrization 0, 
l22j in the resonance region (W < 2.5 GeV), and a QCD- 
like fit from Ref. H3 in the DIS (IF > 2.5 GeV): 



F 2 (x,Q 



F 2 (x,Q 2 ) 
F 2 B (x,Q 2 ) 



F 2 B {x,Ql) 



(25) 
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FIG 

^,2 



7: Integrands of the inelastic Nachtmann moments at 
0.825 GeV 2 : circles represent the integrand of the Mi\ 
squares show the integrand of the M4; triangles - the inte- 
grand of the Ma; crosses - the integrand of the Ms- 



The difference between the real and bin-centered data 



rcent 
°F 2 



(x,Q') = F 2 (x,Q') 



1 



F B (x,Q 2 ) 



(26) 



was added to the systematic errors of F 2 in extracting 
the Nachtmann moments. As an example, Fig. [7\ shows 
the integrands of the first four moments as a function of 
x at fixed Q 2 . The significance of the large x region for 
various moments can clearly be seen . 

To have a data set dense in x, which reduces the error 
in the numerical integration, we performed an interpola- 
tion, at each fixed Qq, when two contiguous experimen- 
tal data points differed by more than V. The value of V 
depended on kinematics; in the resonance region where 
the structure function exhibits strong variations, V had 
to be smaller than half the resonance widths and was 
parametrized as V = 0.04/ 1 + i/Q^ld . Ab ove the 

resonances, where F 2 is smooth, we only accounted for 
the fact that the available x region decreases with de- 
creasing Q 2 (V = 0.1 1 + y/Q 2 /lO ). Finally in the low 

x region (x < 0.03) where the F 2 shape depends weakly 
on Q 2 , but strongly on x, we set V = 0.015. Chang- 
ing these V values by as much as a factor two produced 
changes in the moments that were much smaller than the 
systematic errors. 

To fill the gap within two contiguous points x a and Xb, 
we used the interpolation function F 2 nt (x, Qq) defined as 
the parametrization from Ref. normalized to the ex- 
perimental data on both edges of the interpolating range. 
Assuming that the shape of the fit is correct: 

F 2 int (x,Q 2 ) = a(Ql)F 2 B (x,Ql) , (27) 



where the normalization factor a(Qo) is defined as the 
weighted average, evaluated using all experimental points 



«(Qo) = &(Qi>) 

\xj — x b I <A 



\xi —x a I < A 



E 



F 2 (x t ,Q 2 )/F 2 B ( Xl ,Q 2 ) 
(6% at ( Xl ,Q 2 )) 2 



E 

i 



F 2 (x J: Q 2 )/F B { Xjl Q 2 ) 
(Spfix^Q 2 )) 2 



(28) 



where 5p^ lt (xj, Qq) is the statistical error relative to F 2 B 



and 



s N (Ql) = 



\xi-x a \<A 

E 



(«Sf°*(*i.Q§))' 



1 



\xj— xf,|<A 

V - 



-1/2 



(29) 



is the statistical uncertainty of the normalization. There- 
fore, the statistical error of the moments calculated ac- 
cording the trapezoidal rule 0] was increased by adding 
the linearly correlated contribution from each interpola- 
tion interval as follows: 

x b 

f fn+1 

6™ rm (Ql)= S N (Q 2 ) / dx^Fi^Q 2 ) 



3 + 3(n + l)r + n(n + 2)r 2 

(n + 2)(n + 3) ' ( ' 

Since we average the ratio F 2 (xi, Q 2 ,)/F 2 B (xi, Qq), A 
is not affected by the resonance structures, and its value 
was fixed to have more than two experimental points in 
most cases; therefore, A was chosen equal to 0.03 in the 
resonance and in the very low x regions, and to 0.05 in 
the DIS region. In Fig. [S] we show how this interpolation 
is applied: the thin continuous line represents the original 
function F 2 B (x, Q 2 ) and the heavy continuous line repre- 
sents the result of the interpolation F 2 nt (x, Q 2 ). We also 
checked that the moments do not show any dependence 
on the A values. 

To fill the gap between the last experimental point 
and one of the integration limits (x a — or ij = 1) 
we performed an extrapolation at each fixed Qq us- 
ing F B (x, Q 2 ) including its systematic error given in 
Ref. [l|. 

As an extension of the analysis, the world data at Q 2 
above 5 (GeV/c) 2 were analyzed in the same way as de- 
scribed above. The only differences were the Q 2 bin size, 
which was chosen equal to 5% of Q 2 , and the values of 
the parameters V and A. In addition, the bins were sit- 
uated not continuously, but only where data exist. Since 
at large Q 2 the shape of F 2 (x) remains almost constant 
with changing Q 2 both parameters V and A were fixed: 
in the resonance region (W < 1.8 GeV) to value 0.01; 
0.1 in the DIS; and V = 0.005 and A = 0.01 at very low 
x (x < 0.03). The results for AP n n (Q 2 ) did not exhibit 
any significant dependence on the choice of the parame- 
ter values. The results are reported together with their 
statistical and systematic errors in Table ITT1 
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FIG. 8: Example of the interpolating procedure. The mean- 
ing of the curves and symbols is described in the text. 



H. Systematic Errors of the Moments 

The systematic error consists of genuine uncertainties 
in the data given in Refs. S El J2§ M EJ HI IH 
SIllllllSlliElilElilililand uncertain- 
ties in the evaluation procedure. To estimate the first 
type of error we had to account for using many data sets 
measured, in different laboratories, and with different de- 
tectors. In the present analysis we assume that different 
experiments are independent and therefore only system- 
atic errors within one data set are correlated. 

Thus, an upper limit for the contribution of the sys- 
tematic error from each data set was evaluated in the 
following way: 

• we first applied a simultaneous shift to all exper- 
imental points in this set by an amount equal to 
their systematic error; 

• then the inelastic n-th moments obtained using 
these distorted data M^^(Q 2 ) were compared to 

the original moments M™(Q 2 ) evaluated with no 
systematic shifts; 

• finally the deviations for each data set were 
summed in quadrature as independent values: 



M 



l n (Q 2 )\ 



N S / 



(31) 

where Ns is the number of available data sets. 
The resulting error was summed in quadrature to 
5™ rm {Q 2 ) to finally evaluate the total systematic 
error of the n-th moment. 

The second type of error is related to the bin center- 
ing, interpolation and extrapolation. The bin centering 



systematic uncertainty was estimated as follows: 

$n(Q 2 ) =Y,K n (x i ,tf)w i (0*)8%*(x il Cf) , (32) 



where according the Nachtmann moment definition and 
the trapezoidal integration rule: 



K n {x l7 Q 2 ) 



3 + 3(n + l)n + n(n + 2)r| 
xf (n + 2)(n + 3) 



WiiQ 2 ) = (x i+ i - a?i_i)/2 . 



(33) 



The relative systematic error of the interpolation was 
estimated as the possible change of the fitting function 
slope in the interpolation interval, and it was evaluated 
as a difference in the normalization at different edges: 



5s(Ql) = 



1 N ^ <A F 2 { Xi ,Ql) 
Ni F?{x u Ql) 



N 



1 '^< A F 2 ( Xj ,Ql) 



F 2 (xj, Qq) 



(34) 



where Ni and Nj are the number of points used to eval- 
uate the sums. Since the structure function F 2 (x,Q 2 ) 
is a very smooth function of x below resonances, on the 
limited x-interval (smaller than V) the linear approxi- 
mation gives a good estimate. Thus, the error given in 
Eq. El accounts for such a linear mismatch between the 
fitting function and the data on the interpolation inter- 
val. Meanwhile, the CLAS data cover all the resonance 
region and no interpolation was used there. Therefore, 
the total systematic error introduced in the correspond- 
ing moment by the interpolation can be estimated as 



#(Qo)= S s (Qo) / dx^—F 2 B {x,Ql) 



i n+1 

x° 



3 + 3(n + l)r + n(n + 2)r 2 
(n + 2)(n + 3) 



(35) 



The systematic errors obtained by these procedures 
were summed in quadrature: 



€{Q 2 ) = ^K(Q2 W + {S C m) 2 + {S I im) 2 . (36) 

In order to study the systematic error on the extrapola- 
tion at very low x we have performed a test of the func- 
tional form dependence comparing moments presented 
here with those obtained by using the fitting function 
from the neural network parametrization of Ref. |52| . 
The difference is significant only for M 2 and it appeared 
to be smaller than 5^(Q 2 ) given by Eq. ESI (see Fig.EJ). 
The difference was added to S^(Q 2 ) in quadrature to 
evaluate the total systematic error of the n-th moment. 
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FIG. 9: Errors of the inelastic Nachtmann moment M2 in 
percentage: empty circles represent statistical errors; empty 
crosses show the systematic error obtained in Eq. 1361 the dif- 
ference between inelastic moments extracted using two dif- 
ferent F2 parametrizations fl4| and [H^l at W > 2.5 GeV is 
shown by stars. 
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Q~ (GeV/c) 

FIG. 10: The inelastic Nachtman moments extracted from 
the world data, including the new CLAS results, are shown 
as the solid shapes ; while the solid line represents moments 
obtained in Ref. flijl. The empty crosses indicate the Nacht- 
man moments determined in Ref. Q. Errors are statistical 
only. 



IV. EXTRACTION OF LEADING AND 
HIGHER TWISTS 



In this Section, we present our twist analysis of the mo- 
ments we have extracted, which are presented in Fig. 1101 
As already shown in Refs. 01 an d the extraction 
of higher twists at large x depends significantly on the 
effects of pQCD high-order corrections. In particular, 
the use of the well established NLO approximation for 
the leading twist is known to lead to unreliable results 
for the higher twists [l^. Therefore, hereafter we fol- 
low Ref. [lij, where the pQCD corrections beyond the 
NLO are estimated according to soft gluon resummation 
(SGR) techniques. 

As far as power corrections are concerned, sev- 
eral higher-twist operators exist and mix under the 
renormalization-group equations. Such mixings are 
rather involved and the number of mixing operators in- 
creases with the order n of the moment. Since a com- 
plete calculation of the higher-twist anomalous dimen- 
sions is not yet available, we use the same phenomenolog- 
ical ansatz already adopted in Refs. and 0|. Thus, 
our extracted Nachtmann moments are analyzed in terms 
of the following twist expansion 



M»(Q^= Vn (Q')+HT n (Q*) , 



(37) 



where r] n (Q 2 ) is the leading twist moment and HT n (Q 2 ) 
is the higher-twist contribution given by (53j 



HT n (Q 2 ) 




(38) 



where the logarithmic pQCD evolution of the 
twist-r contribution is accounted for by the term 
[a s (Q 2 ) / asip 2 )] 1 ^ (corresponding to the Wilson coeffi- 
cient E nT {n,Q 2 ) in Eq. EJ) with an effective anomalous 
dimension and the parameter (equal to the 
matrix element O nT (/i) in Eq. |3J) represents the overall 
strength of the twist-r term at the renormalization scale 
Q 2 = ii 2 . In Eqs. 1371381 four higher- twist parameters 
appear, while the twist-2 moment i] n {Q 2 ) is generally 
given by the sum of a non-singlet and singlet terms, 
leading to three unknown parameters, namely the values 
of the gluon, non-singlet and singlet quark moments 



at the factorization scale Q 



However, the 



decoupling in the pQCD evolution of the singlet quark 
and gluon densities at large x allows one to consider 
a pure non-singlet evolution for n > 4 (cf. This 
means that we have only one twist-2 parameter for 
n > 4, namely the value of the twist-2 moment rj n (fj, 2 ) 
at the factorization scale Q 2 = /i 2 . The resummation of 
soft gluons does not introduce any further parameter in 
the description of the leading twist. Explicitly, for n > 4 
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the leading twist moment rj n (Q 2 ) is given by 



Vn{Q 2 ) = A 



a s (Q 2 ) 



a s (fi 2 
/ I", , as{Q 2 ) ^(NLO) 

R; 



,G„(Q 2 )_ 



Ots(Q 2 ) P NS 



4tt 



(39) 



where the quantities 7^ s , C^jg°^ and R^ s can be read 
off from Ref. [T^j. In Eq.Eithe function G n (Q 2 ) is the 
key quantity of the soft gluon resummation. At next-to- 
leading-log it reads as 

G n (Q 2 ) = ln(n)Gi(A„) + G 2 (A„) + O^ln^H] ,(40) 

where A„ = Poa s (Q 2 )ln(n)/4TT and 



Gi(A) 
G 2 (A) 



G F — [A + (1-A)ln(l-A)] 

-G F — ln(l - A) 

Po 



-G F ^[A + ln(l-A)] 



8_K 

4/?i 
Pi 



A + ln(l - A) + -hi (1- A) 



(41) 



with C F = (N 2 - l)/(2iV c ), k = A c (67/18 - tt 2 /6) - 
5iV//9, /3 = 11 - 2JV//3, and iV> being the number of 
active flavors. Note that the function G 2 (A) is divergent 
for A — ► 1; this means that at large n (i.e. large x) the 
soft gluon resummation cannot be extended to arbitrarily 
low values of Q 2 . Therefore, for a safe use of present SGR 
techniques we will work far from the above-mentioned di- 
vergences by limiting our analyses of low-order moments 
(n < 8) to Q 2 > 1 (GeV/c) 2 

All the unknown parameters, namely the twist-2 
parameter A n as well as the higher-twist parameters 
On 1 , In^ 1 , &n ; , were simultaneously determined from 
a ^-minimization procedure in the Q 2 range between 1 
and 100 (GeV/c) 2 . In such a procedure only the sta- 
tistical errors of the experimental moments were taken 
into account, as well as the updated Particle Data Group 
value a s (M§) = 0.118 jH^], and a renormalization scale 
equal to [i 2 = 10 (GeV/c) 2 . The uncertainties of the var- 
ious twist parameters were then obtained by adding the 
systematic errors to the experimental moments and by 
repeating the twist extraction procedure. The parame- 
ter values are reported in Table ITTT1 where it can be seen 
that: the leading twist is determined with a few percent 
uncertainty, while the precision of the extracted higher 
twists increases with n reaching an overall 10% for n = 6 
and 8, thanks to the remarkable quality of the CLAS 
data at large x. Note that the leading twist is directly 
extracted from the data, which means that no specific 
functional shape of the parton distributions is assumed 



in our analysis. The contribution of higher twists to M 2 
was too small to be extracted by the present procedure. 

Our results, including the uncertainties for each twist 
term separately, are reported in Fig. II II for n > 2, while 
the ratio of the total higher-twist contribution to the 
leading twist is shown in Fig. ^| In addition, the ex- 
tracted leading twist contribution is reported in Table llVl 
It can be seen that: 

• the extracted twist-2 term yields an important con- 
tribution in the whole Q 2 -range of the present anal- 
ysis; 

• the Q 2 -behaviour of the data leaves room for a 
higher-twist contribution positive at large Q 2 and 
negative at Q 2 ~ 1 — 2 (GeV/c) 2 ; such a change of 
sign requires in Eq. |2H1 a twist-6 term with a sign 
opposite to that of the twist-4. As already noted 
in Refs. 0, ^| , such opposite signs make the total 
higher-twist contribution smaller than its individ- 
ual terms; 

• the total higher-twist contribution is significant at 
Q 2 w few (GeV/c) 2 , but it is less than ~ 20% of 
the leading twist for Q 2 > 5 (GeV/c) 2 . 



V. CONCLUSIONS 

We extracted the Fi structure function in a contin- 
uous two-dimensional range of Q 2 and x from the in- 
clusive cross section measured with the CLAS detector. 
Using these data, together with the previously available 
world data set, we evaluated the Nachtmann moments 



M 2 {Q 2 



M 4 (Q 2 ,x), M 6 (Q 2 ,x) and M 8 (Q 2 ,x) in the 



Q 2 range 0.05 - 100 (GeV/c) 2 . The present data set cov- 
ers a large interval in x thus reducing the uncertainties in 
the integration procedure. The Nachtmann moments ob- 
tained in this work have been analyzed in terms of a twist 
expansion in order to simultaneously extract both the 
leading and the higher twists. The former has not been 
treated at fixed order in perturbation theory, but higher- 
order corrections of pQCD were taken into account by 
means of soft gluon resummation techniques. Higher 
twists have been treated phenomenologically by intro- 
ducing effective anomalous dimensions. The range of the 
analysis was quite large, ranging from 1 to 100 (GeV/c) 2 . 
The leading twist is determined with a few percent uncer- 
tainty, while the precision of the higher twists increases 
with n reaching an overall 10% for n = 6 and 8, thanks 
to the remarkable quality of the experimental moments. 

The main results of our twist analysis can be summa- 
rized as follows: i) the contribution of the leading twist 
calculated in the frame of pQCD at NLO remains dom- 
inant down to 2Q 2 /n ~ 1 (GeV/c) 2 , where n is the mo- 
ment order. This leads to the conclusion that a pQCD- 
based description of the proton structure is relevant also 
at low Q 2 , with significant but not crucial corrections; 
ii) the total contribution of the multi-parton correlation 
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TABLE II: The inelastic Nachtmann moments for n = 2, 4, 6 and 8 evaluated in the interval 0.05 < Q 2 < 100 (GeV/c) 2 . The 
moments are labeled with * when the contribution to the integral by the experimental data is between 50% and 70%; all other 
values were evaluated with more than 70% data coverage. The data are reported together with the statistical and systematic 
errors. 



m- ,<;,\ c y 


M 2 (Q")xl0- i 


M 4 (Q^)xl0^ 


M 6 (Q>10- J 




0.075 


0.202 ± 0.002 ± 0.009 


0.016 ± 0.0005 ± 0.001 


0.0019 ± 0.0001 ± 0.0001 




0.125 


0.451 ± 0.006 ± 0.025 


0.072 ± 0.002 ± 0.004 


0.017 ± 0.001 ± 0.001 




0.175 


0.638 ± 0.005 ± 0.025 


0.162 ± 0.002 ± 0.007 


0.060 ± 0.001 ± 0.003 


0.0025 ± 0.0001 ± 0.0001 


0.225 


0.775 ± 0.003 ± 0.026 


0.248 ± 0.001 ± 0.008 


0.119 ± 0.001 ± 0.004 


0.0064 ± 0.0001 ± 0.0002 


0.275 


0.910 ± 0.004 ± 0.030 


0.364 ± 0.002 ± 0.015 


0.218 ± 0.002 ± 0.010 


0.0146 ± 0.0001 ± 0.0007 


0.325 


1.000 ± 0.002 ± 0.040 


0.465 ± 0.0005 ± 0.026 


0.328 ± 0.0005 ± 0.020 


0.0259 ± 0.00005 ± 0.0017 


0.375 


1.114 ± 0.002 ± 0.047 


0.587 ± 0.0005 ± 0.033 


0.478 ± 0.0005 ± 0.030 


0.0439 ± 0.00005 ± 0.0029 


0.425 


1.209 ± 0.005 ± 0.037 


0.704 ± 0.001 ± 0.034 


0.644 ± 0.001 ± 0.038 


0.0670 ± 0.0001 ± 0.0043 


0.475 


1.298 ± 0.008 ± 0.036 


0.839 ± 0.003 ± 0.023 


0.858 ± 0.003 ± 0.024 


0.1002 ± 0.0004 ± 0.0030 


0.525 


1.347 ± 0.004 ± 0.047 


0.916 ± 0.003 ± 0.038 


1.010 ± 0.005 ± 0.046 


0.1279 ± 0.0008 ± 0.0062 


0.575 


1.419 ± 0.003 ± 0.049 


1.023 ± 0.002 ± 0.050 


1.215 ± 0.002 ± 0.068 


0.1660 ± 0.0003 ± 0.0101 


0.625 


1.444 ± 0.006 ± 0.059 


1.110 ± 0.003 ± 0.041 


1.413 ± 0.005 ± 0.057 


0.2079 ± 0.0009 ± 0.0090 


0.675 


1.514 ± 0.004 ± 0.051 


1.191 ± 0.001 ± 0.062 


1.603 ± 0.002 ± 0.098 


0.2507 ± 0.0005 ± 0.0168 


0.725 


1.554 ± 0.006 ± 0.050 


1.267 ± 0.001 ± 0.059 


1.785 ± 0.002 ± 0.102 


0.2946 ± 0.0004 ± 0.0190 


0.775 


1.578 ± 0.007 ± 0.049 


1.345 ± 0.002 ± 0.053 


1.996 ± 0.002 ± 0.087 


0.3484 ± 0.0005 ± 0.0160 


0.825 


1.606 ± 0.006 ± 0.050 


1.389 ± 0.002 ± 0.066 


2.130 ± 0.003 ± 0.117 


0.3860 ± 0.0006 ± 0.0233 


0.875 


1.625 ± 0.019 ± 0.074 


1.452 ± 0.005 ± 0.065 


2.320 ± 0.004 ± 0.122 


0.4393 ± 0.0008 ± 0.0254 


0.925 


1.649 ± 0.014 ± 0.040 


1.500 ± 0.005 ± 0.058 


2.476 ± 0.005 ± 0.119 


0.4866 ± 0.0010 ± 0.0264 


0.975 


1.669 ± 0.013 ± 0.044 


1.553 ± 0.005 ± 0.058 


2.651 ± 0.007 ± 0.113 


0.5416 ± 0.0015 ± 0.0254 


1.025 


1.673 ± 0.011 ± 0.049 


1.584 ± 0.004 ± 0.061 


2.785 ± 0.011 ± 0.116 


0.5887 ± 0.0030 ± 0.0248 


1.075 


1.706 ± 0.011 ± 0.046 


1.597 ± 0.004 ± 0.067 


2.820 ± 0.005 ± 0.140 


0.6048 ± 0.0012 ± 0.0322 


1.125 




1.648 ± 0.003 ± 0.076 


3.002 ± 0.005 ± 0.150 


0.6627 ± 0.0013 ± 0.0370 


1.175 




1.701 ± 0.004 ± 0.055 


3.179 ± 0.007 ± 0.117 


0.7236 ± 0.0018 ± 0.0298 


1.225 


1.722 ± 0.009 ± 0.045 


1.706 ± 0.005 ± 0.066 


3.245 ± 0.009 ± 0.154 


0.7525 ± 0.0020 ± 0.0402 


1.275 


1.736 ± 0.006 ± 0.086 


1.732 ± 0.005 ± 0.060 


3.364 ± 0.012 ± 0.126 


0.8021 ± 0.0036 ± 0.0309 


1.325 


1.792 ± 0.015 ± 0.050 


1.828 ± 0.004 ± 0.076 


3.561 ± 0.011 ± 0.178 


0.8556 ± 0.0033 ± 0.0475 


1.375 


1.798 ± 0.027 ± 0.055 


1.839 ± 0.004 ± 0.082 


3.630 ± 0.008 ± 0.189 


0.8864 ± 0.0024 ± 0.0516 


1.425 


1.815 ± 0.007 ± 0.049 


1.873 ± 0.004 ± 0.073 


3.741 ± 0.011 ± 0.173 


0.9280 ± 0.0032 ± 0.0492 


1.475 


1.833 ± 0.006 ± 0.053 


1.899 ± 0.004 ± 0.073 


3.839 ± 0.010 ± 0.154 


0.9669 ± 0.0031 ± 0.0397 


1.525 


1.844 ± 0.008 ± 0.055 


1.931 ± 0.004 ± 0.082 


3.968 ± 0.012 ± 0.187 


1.0158 ± 0.0042 ± 0.0488 


1.575 


1.833 ± 0.006 ± 0.065 


1.940 ± 0.004 ± 0.096 


4.022 ± 0.010 ± 0.249 


1.0395 ± 0.0033 ± 0.0725 


1.625 


1.862 ± 0.020 ± 0.053 


1.953 ± 0.005 ± 0.091 


4.116 ± 0.010 ± 0.252 


1.0859 ± 0.0034 ± 0.0772 


1.675 




1.957 ± 0.005 ± 0.083 


4.170 ± 0.011 ± 0.231 


1.1173 ± 0.0036 ± 0.0740 


1.725 


1.857 ± 0.023 ± 0.049 


1.998 ± 0.005 ± 0.075 


4.316 ± 0.013 ± 0.218 


1.1680 ± 0.0043 ± 0.0726 


1.775 


1.884 ± 0.063 ± 0.054 


2.020 ± 0.011 ± 0.072 


4.412 ± 0.012 ± 0.194 


1.2081 ± 0.0043 ± 0.0628 


1.825 


1.862 ± 0.010 ± 0.053 


2.024 ± 0.006 ± 0.072 


4.459 ± 0.015 ± 0.168 


1.2338 ± 0.0050 ± 0.0462 


1.875 


1.837 ± 0.015 ± 0.060 


2.014 ± 0.006 ± 0.101 


4.446 ± 0.015 ± 0.256 


1.2363 ± 0.0046 ± 0.0798 


1.925 




2.026 ± 0.006 ± 0.093 


4.551 ± 0.015 ± 0.243 


1.2903 ± 0.0047 ± 0.0755 


1.975 


1.866 ± 0.010 ± 0.059 


2.027 ± 0.007 ± 0.091 


4.539 ± 0.018 ± 0.253 


1.2857 ± 0.0058 ± 0.0824 


2.025 


1.831 ± 0.014 ± 0.046 


2.037 ± 0.007 ± 0.092 


4.677 ± 0.020 ± 0.263 


1.3480 ± 0.0069 ± 0.0867 


2.075 




2.046 ± 0.008 ± 0.084 


4.699 ± 0.022 ± 0.232 


1.3694 ± 0.0084 ± 0.0750 


2.125 


1.870 ± 0.022 ± 0.052 


2.074 ± 0.008 ± 0.092 


4.825 ± 0.022 ± 0.269 


1.4239 ± 0.0082 ± 0.0903 


2.175 


1.846 ± 0.013 ± 0.059 


2.064 ± 0.010 ± 0.098 


4.850 ± 0.024 ± 0.282 


1.4421 ± 0.0092 ± 0.0945 


2.225 




2.053 ± 0.012 ± 0.089 


4.825 ± 0.024 ± 0.267 


1.4442 ± 0.0093 ± 0.0912 


2.275 


1.852 ± 0.020 ± 0.050 


2.062 ± 0.008 ± 0.095 


4.852 ± 0.023 ± 0.271 


1.4606 ± 0.0092 ± 0.0917 


2.325 


1.859 ± 0.012 ± 0.058 


2.081 ± 0.009 ± 0.108 


4.984 ± 0.025 ± 0.291 


1.5149 ± 0.0098 ± 0.0959 


2.375 


1.867 ± 0.012 ± 0.055 


2.060 ± 0.008 ± 0.101 


4.876 ± 0.023 ± 0.275 


1.4832 ± 0.0091 ± 0.0921 


2.425 




2.051 ± 0.008 ± 0.107 


4.860 ± 0.023 ± 0.257 


1.4835 ± 0.0089 ± 0.0850 


2.475 


1.793 ± 0.068 ± 0.089 


2.056 ± 0.010 ± 0.082 


4.971 ± 0.020 ± 0.234 


1.5362 ± 0.0079 ± 0.0796 


2.525 


1.845 ± 0.031 ± 0.066 


2.035 ± 0.010 ± 0.110 


4.899 ± 0.018 ± 0.260 


1.5176 ± 0.0063 ± 0.0751 


2.575 


1.841 ± 0.019 ± 0.052 


2.050 ± 0.010 ± 0.103 


4.972 ± 0.021 ± 0.280 


1.5556 ± 0.0078 ± 0.0915 


2.625 




2.035 ± 0.011 ± 0.122 


4.884 ± 0.024 ± 0.293 


1.5218 ± 0.0087 ± 0.0933 


2.675 




2.018 ± 0.009 ± 0.024 


4.896 ± 0.022 ± 0.277 


1.5457 ± 0.0091 ± 0.0988 


2.725 




2.028 ± 0.011 ± 0.099 


4.933 ± 0.025 ± 0.283 


1.5634 ± 0.0094 ± 0.0970 


2.775 




2.028 ± 0.017 ± 0.107 


4.931 ± 0.029 ± 0.293 


1.5677 ± 0.0096 ± 0.0989 


2.825 


1.836 ± 0.026 ± 0.061 


2.031 ± 0.014 ± 0.118 


5.004 ± 0.028 ± 0.326 


1.6030 ± 0.0098 ± 0.1081 


2.875 


1.839 ± 0.016 ± 0.057 


2.019 ± 0.013 ± 0.108 


4.976 ± 0.027 ± 0.309 


1.6032 ± 0.0098 ± 0.1036 
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Q" [(GoV/c) z j 


M 2 (<3 z )xl0 _1 


M^Q^xlO^ 


M 6 (Q A )xlQ-- > 


M 8 (Q")xlO"' i 


2.925 




2.023 ± 0.016 ± 0.112 


5.007 ± 0.033 ± 0.303 


1.6219 ± 0.0104 ± 0.0970 


2.975 


1.843 ± 0.033 ± 0.050 


2.018 ± 0.014 ± 0.102 


4.983 ± 0.027 ± 0.294 


1.6145 ± 0.0092 ± 0.0941 


3.025 


1.816 ± 0.068 ± 0.058 


1.978 ± 0.016 ± 0.104 


4.926 ± 0.027 ± 0.314 


1.6042 ± 0.0090 ± 0.1057 


3.075 


1.804 ± 0.060 ± 0.055 


1.992 ± 0.022 ± 0.114 


4.942 ± 0.040 ± 0.352 


1.6136 ± 0.0119 ± 0.1222 


3.125 








1.6293 ± 0.0118 ± 0.1491 


3.175 




2.011 ± 0.031 ± 0.141 


5.002 ± 0.064 ± 0.372 


1.6524 ± 0.0159 ± 0.1310 


3.225 




1.968 ± 0.021 ± 0.112 


4.916 ± 0.040 ± 0.358 


1.6289 ± 0.0122 ± 0.1315 


3.275 




2.007 ± 0.022 ± 0.116 


4.985 ± 0.043 ± 0.351 


1.6478 ± 0.0127 ± 0.1304 


3.325 


1.808 ± 0.033 ± 0.080 


1.979 ± 0.014 ± 0.096 


4.944 ± 0.032 ± 0.332 


1.6321 ± 0.0120 ± 0.1264 


3.375 


1.804 ± 0.029 ± 0.055 


1.981 ± 0.016 ± 0.086 


4.976 ± 0.031 ± 0.312 


1.6543 ± 0.0119 ± 0.1243 


3.425 






5.034 ± 0.035 ± 0.325 


1.6787 ± 0.0125 ± 0.1285 


3.475 




1.943 ± 0.013 ± 0.064 


4.915 ± 0.032 ± 0.253 


1.6489 ± 0.0118 ± 0.1079 


3.525 




1.951 ± 0.021 ± 0.088 


4.999 ± 0.052 ± 0.316 


1.6918 ± 0.0168 ± 0.1238 


3.575 






5.021 ± 0.043 ± 0.268 


1.6858 ± 0.0142 ± 0.1049 


3.625 










3.675 




1.930 ± 0.024 ± 0.040 


4.857 ± 0.063 ± 0.310 


1.6493 ± 0.0197 ± 0.1199 


3.725 








1.6698 ± 0.0160 ± 0.1276 


5.967 




1.810 ± 0.015 ± 0.116 


4.597 ± 0.044 ± 0.553 




7.268 




1.743 ± 0.011 ± 0.044 






7.645 






4.374 ± 0.044 ± 0.098 


1.5659 ± 0.0202 ± 0.0396 


8.027 






4.279 ± 0.027 ± 0.135 


1.5205 ± 0.0107 ± 0.0642 


8.434 




1.653 ± 0.014 ± 0.084 


4.223 ± 0.032 ± 0.109 


1.5264 ± 0.0122 ± 0.0419 


8.857 


1.723 ± 0.015 ± 0.041 


1.645 ± 0.019 ± 0.027 


4.108 ± 0.042 ± 0.109 


1.4712 ± 0.0138 ± 0.0566 


9.781 




1.653 ± 0.010 ± 0.061 


4.130 ± 0.034 ± 0.146 


1.4818 ± 0.0167 ± 0.0666 


10.267 


1.752 ± 0.015 ± 0.052 


1.622 ± 0.019 ± 0.031 


4.016 ± 0.035 ± 0.095 


1.4321 ± 0.0113 ± 0.0367 


10.793 






3.987 ± 0.106 ± 0.761 


1.4256 ± 0.0175 ± 0.1103 


11.345 


1.731 ± 0.016 ± 0.041 


1.573 ± 0.018 ± 0.035 


3.853 ± 0.041 ± 0.140 


1.3644 ± 0.0176 ± 0.0793 


11.939 


1.759 ± 0.008 ± 0.042 


1.596 ± 0.013 ± 0.031 


3.910 ± 0.040 ± 0.111 


1.3860 ± 0.0181 ± 0.0574 


13.185 




1.525 ± 0.016 ± 0.021 


3.681 ± 0.029 ± 0.067 


1.3011 ± 0.0091 ± 0.0336 


15.310 


1.686 ± 0.014 ± 0.074 


1.471 ± 0.019 ± 0.032 


3.533 ± 0.044 ± 0.133 




16.902 


1.718 ± 0.010 ± 0.051 


1.450 ± 0.017 ± 0.025 


3.392 ± 0.058 ± 0.073 


1.1752 ± 0.0252 ± 0.0283 


18.697 


1.679 ± 0.033 ± 0.097 


1.401 ± 0.013 ± 0.027 


3.275 ± 0.039 ± 0.088 


1.1377 ± 0.0147 ± 0.0346 


19.629 








*1.1061 ± 0.0221 ± 0.0473 


21.625 


1.677 ± 0.008 ± 0.041 








24.192 


*1.711 ± 0.007 ± 0.114 


1.385 ± 0.008 ± 0.024 




*1.0751 ± 0.0143 ± 0.0433 


26.599 


*1.665 ± 0.062 ± 0.135 








28.192 


*1.702 ± 0.009 ± 0.140 


1.344 ± 0.007 ± 0.037 




*1.0109 ± 0.0096 ± 0.0808 


31.858 


"T.703 ± 0.010 ± 0.096 








36.750 


*1.696 ± 0.013 ± 0.111 


1.314 ± 0.009 ± 0.057 


2.971 ± 0.027 ± 0.313 


*1.0027 ± 0.0135 ± 0.1906 


44.000 


*1.681 ± 0.013 ± 0.085 








49.750 


*1.658 ± 0.019 ± 0.101 








57.000 


*1.694 ± 0.017 ± 0.170 








64.884 


*1.636 ± 0.043 ± 0.114 


1.222 ± 0.053 ± 0.044 


2.708 ± 0.082 ± 0.193 


*0.8945 ± 0.0161 ± 0.1164 


75.000 




*1.206 ± 0.008 ± 0.025 


*2.651 ± 0.024 ± 0.150 




88.000 


*1.669 ± 0.088 ± 0.075 


*1.199 ± 0.038 ± 0.035 


*2.630 ± 0.057 ± 0.202 




99.000 




"T.179 ± 0.012 ± 0.034 


*2.568 ± 0.029 ± 0.228 





TABLE III: Extracted parameters of the twist expansion. The contribution of higher twists to M 2 was too small to be extracted 
by the present procedure. 





M 2 


Ma 


M 6 


Ms 


A 2 


0.174± 0.006 


(1.61± 0.04)xl0~ 2 


(3.98± 0.18)xl0~ :l 


(1.39± 0.07) xl0" a 


a< 4 > 


(1.4 ± 1.8) xl0~ 3 


(3.6± 1.4) xl0~ 4 


(1.9± 0.14)xl0~ 4 


(1.69± 0.16) xl0~ 4 


7 (4) 




5.7 ± 0.6 


7.4 ± 0.3 


6.2 ± 0.3 


a™ 




(-9.5 ± 3.4)xl0~ 5 


(-6.57 ± 0.53) xl0~ 5 


(-5.75 ± 0.44) xl0~ 5 


7 («) 




4.4 ± 0.6 


5.7 ± 0.3 


4.6 ± 0.3 
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TABLE IV: The extracted leading twist contribution r] n (Q 2 ) (see Ea . I39[l shown in Fig. 1111 reported with systematic errors. 



Q [(GeV/c) J 


T} 2 {Q jxlO 


r]4.(Q JxlO 


„ { \ 1 A - ^ 

ri 6 (Q JxlO 


Vs(Q )xio 


1.025 


o i o I n at 

1.16 ± (J. 07 


o eo i n AA 
0.02 ± 0.U9 


i o c c I n nc 
1.DD0 ± U.U7 


i ooo i f\ nc c 

1.223 ± U.Uoo 


1.0/5 


oil _i_ n nc 
2.11 ± U.Uf 


o a n _i_ n no 
0.49 ± U.U9 


i coo _i_ n nc 
1.022 zh U.U/ 


i noo _i_ n nrr 
1.1)22 ± U.Uoo 


1.125 


o nn i a at 

2.U9 ± U.U7 


o o o i n a o 

0.08 ± U.U8 


1 A 1 n 1 A AA' 

1.41U ± U.Uo 


n o o o i a a ,i 7 

U.88o ± U.U47 


1.1 to 


o no _i_ n nc 
2.U8 ± U.U/ 


o oo _i_ n no 
0.28 ± U.U8 


i o i n _i_ n nc 
l.ol9 ± U.UD 


n C0 1 1 A A A 1 

U. (oi ± U.U41 


1.225 


O A7 1 C\ nC 

2.07 ± 0.07 


OTA 1 A no 

3.19 ± 0.08 


i o a o i a at 

1.24o zb U.Uo 


n cn a i a\ n o t 

0.704 ± 0.037 


1.275 


2.05 ± 0.07 


oil i n no 

3.11 ± 0.08 


1 1 7fl i n A\ C 

1.179 ± 0.05 


n c /t o i n a o i 

0.643 ± 0.034 


1.325 


o a a i n a\ t 

2.04 ± 0.07 


O C\ A i n AT 

3.04 ± 0.07 


i ioc i a n c 

1.125 ± U.U5 


n cno i a noi 

0.593 ± 0.031 


1.375 


o no i a at 

2.Uo ± U.07 


o nc i n nc 

2.97 ± 0.U7 


i A77 i n nc 

1.U77 ± U.U5 


n c c o i f\ n o n 

U.ooo ± U.U29 


1 /I OK 

1.425 


o no _i_ n nc 
2.U2 ± U.U/ 


o m _i_ n nT 
2.91 ± U.U/ 


i noc _i_ n nc 
l.UoD ± U.UO 


n ma i a noc 

u.oiy ± u.U2 ^ 


1.475 


O A1 1 A AC 

2.01 ± 0.07 


o oc i n nT 

2.86 ± 0.07 


n nnn i n n/i 

0.999 ± 0.04 


n a nn i a\ noc 

0.490 ± 0.026 


1.525 


O A A 1 A A"r 

2.00 ± 0.07 


o o i i n nT 

2.81 ± 0.07 


n nAc i n n^i 

0.966 ± 0.04 


n /t a c i n no a 
U.4oo ± U.U24 


1.575 


1 A A 1 A AT 

1.99 ± 0.07 


A T A 1 A AT 

2.76 ± 0.07 


A A O f 1 A A 1 

0.936 ± 0.04 


A A A O 1 A A. "1 ■ > 

0.443 ± 0.023 


1.625 


1 AO 1 A AT 

1.98 ± 0.07 


O CO 1 A\ AC 

2.72 ± 0.07 


n n i n i n n a 

0.910 ± 0.04 


A A O A \ A noo 

U.424 ± U.U22 


1.675 


1 AT 1 A AT 

1.97 ± 0.07 


O A" O 1 A AT 

2.68 ± 0.07 


n o oc i c\ c\ a 

0.886 ± 0.04 


A A AC i c\ no 1 

0.407 ± 0.021 


1.725 


1 CvO 1 A A.T 

1.96 ± 0.07 


O A A 1 A A A' 

2.64 ± 0.06 


A OA/I 1 A A/1 

U.od4 zb U.U4 


n o n o i a\ n o 1 
U.o92 ± U.U21 


1.775 


1 n C 1 A A. t 

1.95 ± 0.07 


O A 1 1 A A A' 

2.61 ± 0.06 


n O A A 1 A A/1 

U.844 zb U.U4 


n o c o i n no n 

0.378 ± 0.020 


1.825 


1 AK I A /"I T 

1.95 ± 0.07 


A P* T 1 A A A" 

2.57 ± 0.06 


A O A C 1 A A 1 

0.825 ± 0.04 


A O / "' / I A A 1 A 

0.366 ± 0.019 


1.875 


1 A A 1 A AT 

1.94 ± 0.07 


O C /I | /"l f"lA 

2.54 ± O.Oo 


n ono i c\ c\ a 

0.808 ± 0.04 


a o c c i n m n 

0.355 ± 0.019 


1.925 


1 AO 1 A A"~ 

1.93 ± 0.07 


O C i 1 A A A' 

2.51 ± O.Oo 


n 7ao i a no 

0.792 ± 0.03 


n O A A 1 A A 1 D 

U.o44 ± U.U18 


1.975 


1 AO 1 A AT 

1.93 ± U.07 


o a n i n n a* 

2.49 ± U.Uo 


r\ i-ri-rrj i A AO 

0.777 ± U.Uo 


nooc i a ai o 
U.OOO ± U.U18 


o no c 


i no _l n nc 
1.92 ± U.Uf 


o jfi _l a nc 
2.4o ± U.Uo 


n ceo _i_ n no 
0. (Do zb U.Uo 


n ooc _l n m c 
U.o2o ± U.U1 


2.075 


i n i i a a t 

1.91 ± 0.07 


O A A 1 n AA' 

2.44 ± U.Uo 


n C C A I A A o 

0.750 ± 0.03 


n o i o i a a i 7 

0.318 ± 0.017 


2.125 


i a i i n at 

1.91 ± 0.07 


O A 1 J n AA' 

2.41 ± U.Uo 


n coo i n no 

0.738 ± 0.03 


A O 1 1 I A A 1 A 

0.311 ± 0.016 


2.175 


1 A A 1 A AT 

1.90 ± 0.07 


o on i n nc 

2.39 ± O.Oo 


n CO A I A AO 

0.726 ± 0.03 


n o n a i n m c 

0.304 ± 0.016 


2.225 


1 A A 1 A AT 

1.90 ± 0.07 


o ot i a nc 

2.37 ± O.Oo 


n CI C 1 A AO 

0.715 ± 0.03 


n ono i n a 1 c 

0.298 ± 0.016 


2.275 


1 OA 1 A AT 

1.89 ± 0.07 


o o cr i a r\ C' 

2.35 ± O.Oo 


A 7AA 1 A AO 

0.706 ± 0.03 


n ono i n ai 1 c 
U.292 ± U.Ulo 


2.325 


1 OA 1 A AT 

1.89 ± 0.07 


O O O 1 A A A' 

2.33 ± O.Oo 


n enc i n no 

0.697 ± 0.03 


nooc i a a i r 

0.287 ± 0.015 


2.375 


1 OA 1 A AT 

1.89 ± 0.07 


O OO 1 f\ A A' 

2.32 ± U.Uo 


neon i a n o 

0.689 ± 0.03 


n o o o i a a i r 

0.283 ± 0.015 


2.425 


1 O O 1 A A "~ 

1.88 ± 0.07 


A OA 1 A\ A A" 

2.30 ± 0.06 


A A' OA 1 A A'"> 

0.682 ± 0.03 


A O T A 1 A\ A 1 i 

0.279 ± 0.014 


2.475 


1 OO 1 A A\t 

1.88 ± 0.07 


o oo i a nc 
2.28 ± U.Uo 


n O C C 1 A A o 

0.675 ± 0.03 


aocc i n m a 

0.275 ± 0.014 


2.525 


1 O O 1 A AC 

1.88 ± 0.06 


o ot i a nc 

2.27 ± 0.06 


n a o o i a a o 

0.668 ± 0.03 


noci i r\ r\~t a 
U.271 ± U.U14 


2.575 


1 OO 1 A C\f 

1.88 ± 0.06 


O O CT 1 a A A' 

2.25 ± U.Uo 


n a c i i n n o 

0.661 ± 0.03 


nocc i r\ r\~\ a 

0.267 ± 0.014 


2.625 


1 07 1 A AA 

1.87 ± 0.06 


O O A 1 A AC 

2.24 ± U.Uo 


n a c /t i a n o 

0.654 ± 0.03 


noc/t i n m A 
U.2o4 ± U.U14 


2.675 


i oc i a nc 

1.87 ± 0.06 


o oo i a n c 
2.23 ± U.Uo 


n a a o i a n o 

0.648 ± 0.03 


nocn i nmo 

0.260 ± 0.013 


2.725 


1 OT I a i \t ' 

1.87 ± 0.06 


o o i i a n c 
2.21 ± U.Uo 


n a /i o i a n o 

0.642 ± 0.03 


nocc i nmo 

0.257 ± 0.013 


2.775 


1 OT 1 A A C 

1.87 ± 0.06 


AAA 1 A AT 

2.20 ± 0.05 


A A O t 1 r\ f\c\ 

0.637 ± 0.03 


A A r 1 1 A A 1 ■"> 

0.254 ± 0.013 


2.825 


1 OC 1 A nc 

1.86 ± 0.06 


O 1 n 1 A AC 

2.19 ± 0.05 


n a o i i n n o 

0.631 ± 0.03 


a o c i i nmo 
U.251 ± U.Ulo 


2.875 


1 OC 1 A A/ - "' 

1.86 ± 0.06 


O 1 O 1 A AC 

2.18 ± 0.05 


n a OA* i a no 

0.626 ± 0.03 


n o a o i nmo 
U.248 zb U.Ulo 


2.925 


i o a i a nc 

1.86 ± 0.06 


o i t i a n c 

2.17 ± 0.05 


n a o i i n n o 

0.621 ± 0.03 


a o a c i nmo 
U.245 ± U.Ulo 


2.975 


i o a i a nc 

1.86 ± 0.06 


o i c i a n c 

2.15 ± 0.05 


n A 1 A 1 A A O 

0.616 ± 0.03 
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3.025 


1.85 ± 0.06 
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O 1 O 1 A A C 

2.12 ± U.Uo 


n c no i n no 

0.602 ± 0.03 


nooc i nmo 

U.2oo ± U.U12 


3.175 


IOC l A A/ ' 

1.85 ± 0.06 


oil i a n c 

2.11 ± U.Uo 


c\ c^rio i a n o 

0.598 ± 0.03 


aooo i nmo 
U.2oo ± U.U12 


3.225 


IOC l a nc 

1.85 ± 0.06 


O 1 A 1 A n c 

2.10 ± 0.05 


n c n o i n n o 

0.593 ± 0.03 


nooi i nmo 
U.2ol ± U.U12 


3.275 


1 O d 1 A A f" 

1.84 ± 0.06 


AAA 1 A AT 

2.09 ± 0.05 


A C O A 1 A A\0 

0.589 ± 0.03 


A AAA 1 A A 1 A 

0.229 ± 0.012 


3.325 


1 O /I 1 A AA' 

1.84 ± 0.06 


O A O 1 A AC 

2.08 ± 0.05 


a r o C 1 A A O 

0.585 ± 0.03 
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0.227 ± 0.012 


3.375 
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O A O 1 A AC 
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O AT 1 A n c 

2.07 ± 0.05 


n CCO 1 A AO 
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1 O A 1 a nc 
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O A C 1 A AC 

2.05 ± 0.05 


n c c i i n n o 

0.571 ± 0.03 


C\ OOO 1 (\ f\ 1 1 

0.220 ± 0.011 


3.575 


1 83 ± 06 


2 04 ± 05 


567 ± 03 


0.218 ± 0.011 


3.625 


1.83 ± 0.06 


2.03 ± 0.05 


0.564 ± 0.02 


0.216 ± 0.011 


3.675 


1.83 ± 0.06 


2.03 ± 0.05 


0.561 ± 0.02 


0.215 ± 0.011 


3.725 


1.83 ± 0.06 


2.02 ± 0.05 


0.558 ± 0.02 


0.213 ± 0.011 


5.967 


1.78 ± 0.06 


1.79 ± 0.04 


0.467 ± 0.02 


0.169 ± 0.009 


7.268 


1.77 ± 0.06 


1.72 ± 0.04 


0.438 ± 0.02 


0.156 ± 0.008 


7.645 


1.76 ± 0.06 


1.70 ± 0.04 


0.431 ± 0.02 


0.153 ± 0.008 


8.027 


1.76 ± 0.06 


1.68 ± 0.04 


0.424 ± 0.02 


0.150 ± 0.008 
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Q z [(GeV/c)'] 








^(Q^xlO^ 


8.434 


1.75 ± 0.06 


1.66 ± 0.04 


0.418 ± 0.02 


0.147 ± 0.008 


8.857 


1.75 ± 0.06 


1.65 ± 0.04 


0.412 ± 0.02 


0.144 ± 0.008 


9.781 


1.74 ± 0.06 


1.61 ± 0.04 


0.400 ± 0.02 


0.139 ± 0.008 


10.267 


1.74 ± 0.06 


1.60 ± 0.04 


0.395 ± 0.02 


0.137 ± 0.007 


10.793 


1.74 ± 0.06 


1.58 ± 0.04 


0.389 ± 0.02 


0.134 ± 0.007 


11.345 


1.73 ± 0.06 


1.57 ± 0.04 


0.384 ± 0.02 


0.132 ± 0.007 


11.939 


1.73 ± 0.06 


1.55 ± 0.04 


0.379 ± 0.02 


0.130 ± 0.007 


13.185 


1.72 ± 0.06 


1.52 ± 0.04 


0.369 ± 0.02 


0.126 ± 0.007 


15.310 


1.71 ± 0.06 


1.48 ± 0.04 


0.355 ± 0.02 


0.120 ± 0.007 


16.902 


1.71 ± 0.06 


1.46 ± 0.04 


0.347 ± 0.02 


0.116 ± 0.006 


18.697 


1.70 ± 0.06 


1.43 ± 0.04 


0.338 ± 0.01 


0.113 ± 0.006 


19.629 


1.70 ± 0.06 


1.42 ± 0.03 


0.334 ± 0.01 


0.111 ± 0.006 


21.625 


1.69 ± 0.06 


1.40 ± 0.03 


0.327 ± 0.01 


0.108 ± 0.006 


24.192 


1.69 ± 0.06 


1.37 ± 0.03 


0.319 ± 0.01 


0.105 ± 0.006 


26.599 


1.68 ± 0.06 


1.35 ± 0.03 


0.312 ± 0.01 


0.102 ± 0.006 


28.192 


1.68 ± 0.06 


1.34 ± 0.03 


0.309 ± 0.01 


0.101 ± 0.006 


31.858 


1.68 ± 0.06 


1.32 ± 0.03 


0.301 ± 0.01 


0.098 ± 0.005 


36.750 


1.68 ± 0.06 


1.29 ± 0.03 


0.293 ± 0.01 


0.095 ± 0.005 


44.000 


1.67 ± 0.06 


1.26 ± 0.03 


0.283 ± 0.01 


0.091 ± 0.005 


49.750 


1.67 ± 0.06 


1.24 ± 0.03 


0.277 ± 0.01 


0.089 ± 0.005 


57.000 


1.67 ± 0.06 


1.22 ± 0.03 


0.270 ± 0.01 


0.086 ± 0.005 


64.884 


1.67 ± 0.06 


1.20 ± 0.03 


0.264 ± 0.01 


0.084 ± 0.004 


75.000 


1.66 ± 0.06 


1.17 ± 0.03 


0.257 ± 0.01 


0.081 ± 0.004 


88.000 


1.66 ± 0.06 


1.15 ± 0.03 


0.250 ± 0.01 


0.078 ± 0.004 


99.000 


1.66 ± 0.06 


1.13 ± 0.03 


0.245 ± 0.01 


0.077 ± 0.004 



effects is not negiigible for Q 2 < 5 (GeV/c) 2 and large 
x corresponding to the resonance region. This can be 
seen by comparing the higher twist contribution to Mg, 
which is more heavily weighted in x, to Afg; iii) dif- 
ferent higher twist terms tend to compensate each other 
in such a way that their sum is small even in a Q 2 re- 
gion where their absolute contributions exceed the lead- 
ing twist. This cancellation is responsible for the duality 
phenomena and leads to the prevailing DIS-inspired pic- 
ture of photon-proton collisions at low Q 2 . 

Therefore, we demonstrated that a precise determina- 
tion of higher twists is feasible with the high quality of 
the new CLAS data. The main limitation of the present 
analysis is the use of a phenomenological ansatz for the 
higher twists. In this respect it is necessary to have bet- 
ter theoretical knowledge of the renormalization group 
behaviour of the relevant higher-twist operators. This 
would directly test QCD in its non-perturbative regime 
through the comparison of predictions obtained from lat- 
tice simulations with these data. 
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APPENDIX A: FIT OF THE RATIO R = a L /a T 

The function R(x, Q 2 ) = <jl/vt was described as: 
R{x,Q 2 ) = 

(l-z) 3 r 0-041gtfc , 0.592 0.331 1 Ti/ ^ o r 

(l-x th y [ C + ~Q 2 ~ ~ (0.09+Q^J W < ' 

(Al) 

0.041g , 0.592 0.331 rovor; 

— h ~TF~ ~ (0.09+Qi) W > / '° 

This parametrization of R(x,Q 2 ) consists of two dif- 
ferent parts: the fit for the DIS region (W > 2.5 
GeV) 0, l5p a nd the function, adjusted to scarce data 
at small in the resonance region (W < 2.5 

GeV). The systematic error on this parametrization was 
estimated according to Ref. as follows: 



0.08 W < 2.5 



0.006C , O01 0.01 wxoc 

— ' ~Cp (0.09+Q 4 ) W > Z -° 



(A2) 



where 



f _i ■ 19 Q 2 0-015625 



l + Q 2 0.015625 + a; 2 
U = aw - 2.5) 
xth = x{W = 2.5) 
All dimensional variables are given in GeV. 
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FIG. 11: Results of the twist analysis. The open squares represent the Nachtman moments obtained in this analysis. The 
solid line is the fit to the moments using Eg. 1371 with the parameters listed in Table UTTI The twist-2, twist-4, twist-6 and higher 
twist (HT) contributions to the fit are indicated. The twist-2 contribution was calculated using Eq. 1391 
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